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ABSTRACT 

The Magellanic Clouds are often considered as outliers in the satellite system 
of the Milky Way because they are irregular and gas-rich galaxies. From their 
large relative motion, they are likely from their first pass near the Milky Way, 
possibly originating from another region of the Local Group or its outskirts. M31 
could have been in a merger stage in its past and we investigate whether or not 
the Large Magellanic Cloud could have been a tidal dwarf expelled during this 
event. Such an hypothesis is tested in the frame of present-day measurements 
and uncertainties of the relative motions of LMC and M31. Our method is to 
trace back the LMC trajectory using several thousands of different configurations 
that sample the corresponding parameter space. 

We find several configurations that let LMC at 50 kpc from M31, 4.3 to 8 
Gyrs ago , depending on the adopted shape of the Milky Way halo. For all 
configurations, the LMC velocity at such a location is invariably slightly larger 
than the escape velocity at such a radius. The preferred solutions correspond to 
a spherical to prolate Milky Way halo, predicting a transversal motion of M31 of 
less than 107 kms _1 and down to values that are close to zero. We conclude that 
from present-day measurements, Magellanic Clouds could well be tidal dwarves 
expelled from a former merger events occurring in M31. 

Subject headings: Galaxies: Local Group - Galaxies: Magellanic Clouds - Galax- 
ies: evolution - Galaxies: dwarf - Galaxies: kinematics and dynamics - Galaxies: 
interactions 
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INTRODUCTION 



The origin of the Magellanic Clouds, as well as the nea r by satellite ga l axies of the Milky 



2009; 


Metz et al. 


2008, 


2009) 



Way (MW), is still a matte r of debate flBesla et al.l 120071 Peebles! 120091 : iKallivavalil et al. 



the Clouds' pro per motions that were car ried out from the Hubble Space Telescope (HST) 



observations by IKallivavalil et al.l ( 12006al Jbl). The total velocity of the Large Magellanic 
Cloud (LMC) in the Galactocentic coordinate is claimed be 378 kms -1 (a transverse velocity 
of ftan = 3 67 kms" 1 an d a ra dial velocity of t> ra d = 89 kms -1 ). Although the revised 
analysis by iPiatek et al.l ( 120081 ) decreases the transverse velocity to 346 kms -1 , both of 



the results 



rom HST data are signifi cantly higher than the previously adopted value, i.e., 



281 kms 1 (Ivan der Marel et al.ll2002l ). At such a high speed, LMC may approach the escape 



velocity at its distance to the MW and the orbital angular momenta is comparable and nearly 
perpendicular to the angular momentum of the MW disk (jKallivayalil et al.ll2009l). This may 
argue for a first passage of the Magellanic Clouds near the Milky Way (see e.g., iBesla et al. 
20071 ). The fact that their morphologies and gas content is at odd with other satellites also 
suggest that they rec ently falls to the Milky Way from the outskirts of the Local Group 
(Ivan den Berghl 120061 ). 



Kallivayalil et al.l (120091 ) construct a model for the Local Group, including Andromeda 
(M31), the MW and the LMC. By solving the equations of motion, they find that M31 may 
have affected the orbit of LMC at a distance of 500-700 kpc about 5 Gyrs ago. Although 
Besla et al.l ( 120071 ) did not investigate the origin of LMC by their models, one interesting 
orbit of LMC is worth to mention here. In their Fig. 14, the orbit under the model of prolate 
MW halo turns close to the direction of M31. The proposition that Magellanic Clouds may 



origin ate fr om M31 was firstly made by iRaychaudhury Lynden-Belll (119891 ). iByrd et al. 
Jl994h and lShuter! Jl992h . 



This Letter revisits this proposition in the frame of the recent discoveries of large 
scale structures surrounding M31 suggesting a very tumultuous past history for this galaxy 
(llbata et al.l 1200 ll . 120041 ; iBrown et al.ll2008l ). We are still lacking of a complete model of 
M31 outskirts although many of its properties are consistent with a past major merger 
(IHammer et al.l 120071 : iBekkil 120101 ) . If true, such a merger sh ould have been gaseous ric h 
enough to allow the reformation of the significant M31 disk (IHammer et al.l 120051 . 120091 ). 
During s uch events gas-rich dwarf g alaxies may be formed from material liberated by the 
collision (jOkazaki fc Taniguchill2000l . and references therein) . It is natural to wonder whether 
or not some tidal dwarf galaxies may have been ejected close to the orbital plane of the hy- 
pothetic merger, which is indeed defined by the actual M31 disk. A significant part of the 
ejected material have angular momentum within small angles from the orbital angular mo- 
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mentum. The M31 disk is seen almost edge-on from the Milky Way, suggesting that the 
Milky Way is located close to the orbital plane of the debris ejected from a major merger 
of M31. This may lead to a fully new interpretation of the Magellanic Clouds, that could 
be tidal dwarves, as massive and concentrated debris lying i n a tidal tail ejected during a 
past event in M31, in the direction of the Milky Way. In fact, lHammer et al.l ( 120101 ) propose 
a major model for the formation of M31 which reproduces most of its properties including 
those of its haunted halo; for some solutions, a significant amount of matter is predicted to 
be ejected from the merger in the direction of the Milky Way. 

The goal of this Letter is to test whether the Magellanic Clouds could have been tidal 
dwarves ejected during a past major merger occurring at the M31 location. The robust 
measurements of LMC prop er motion give a very st rong constrain on its origin by invert- 



ing its past trajectory (e.g., iKallivayalil et al.l 120091 ). On the other hand, there is a large 



uncertainty in the determination of the tangential motion of M31, up to ± 150 km s 1 , (see 



Peebles et al. 2001; Loeb et al. 2005; van der Marel fc Guhathakurta 



20081 1. Taking into ac- 



count all uncertainties, we investigate the possible trajectories of the LMC and whether or 
not it could have approached M31 down to 50 kpc. We solve the equations of motion in 
a dynamical model including MW, M31 and LMC, and throughout the paper, all the 3D 
coordinates, velocities are quoted in the Galactocentric frame that is centered on the MW 
(Ivan der Marel et al.l 120021 ) . We adopt the concordance cosmo logical parameters of H = 70 
km s" 1 Mpc" 1 , Q M = 0.27 and Q A = 0.73. 



2. Analysis 



Fig. [T] shows the 3D positions of the MW, M31 and LMC. A possible unbound trajec- 
tory of LMC in the past is also shown by assuming a zero transverse velocity of M31 for 
our model described below (Sect. 12. ip . This solution is similar to the results presented 
by IKallivayalil et al.l (120091 ) . Given the fact of large uncertainties in the determination 
of tangential motion of M31, up to ±150 kms -1 (see above text), one can expect that 
M31 could have a v x < at present time, meaning that M31 was closer to the past tra- 
jectory of LMC. The 3D velocity of M31 c an be linked to its proper motion on sky by 
following the work by Ivan der Marel et al.l (120021 ) . We adopt the st andard IAU values, 
i.e., Rq = 8.5 kpc and Vq = 220 kms" 1 for the circular velocity (IKerr fc Lynden-Bell 



19861 ). and the solar motion with respect to the local st andard of rest is correcte d by taking 
(£/©, V Ql W q ) = (10.0 ± 0.4, 5.2 ± 0.6, 7.2 ± 0.4) kms" 1 JPehnen fc Binnevlll998h . The basic 
data adopted for M31 are listed in Table [fl as well as the data for LMC. In the following, 
we build a dynamical model of MW, M31 and LMC, then investigate the possible proper 
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Fig. 1. — Upper panels: 3D positions of the MW, M31, LMC. The left panel shows the 
projection in x-y plane while the right panel in y-z plane. The solid lines indicate a possible 
unbound trajectory of LMC in the past. Middle panel: the rot ation curves of M31 . The 
triangles are the measured rotation curve from HI observation (IChemin et al.ll2009l ). The 
red, blue, green and black are the rotation curve of bulge, disk, halo the total of them, 
respectively (Sect. 12.11). Lower panel: the rotation curves of MW. The diamonds are from 



HI observation by lKnapp et al.l (119851 ) . 



Table 1: Basic parameters 
Parameter M31 LMC 
(l,b) (121.174, -21.573) a (280.531, -32.523) b 
Ou w ,m) ( masyr" 1 ) to be investigated (-2.03 ± 0.08, 0.44 ± 0.05) c 
Vfkms" 1 ) -301 d 262.2 b 
D (kpc) 770 e 50.1 f 



a NED. b 


van der Mar el et al. 


( 


2002 


). c ] 


iallivavalil et al. 


( 


2006a 


)■ d 


Courteau h van den Bereh 


(1999). e 


van der Marel & Guhathakurta 


(2008). f 


Kallivavalil et al. 


( 


2009). 



motions of M31 and its impact on the LMC origin. 
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Table 2: Model parameters 



Parameter 


MW 


M31 


LMC 


M virial (1O 12 M ) 


1.0 a 


1.6 a 




^virial (kpc) 


258 b 


300 b 




C c 


15 


18 




M baryon (1O 1O M ) 


5.6 d 


10.9 d 


2.0 


B/T e 


0.15 


0.3 




a b (kpc) 


0.62 


1.0 




r d (kpc) 


2.3 d 


5.8 d 





a Data are from iBesla et al.l (j2007l ) . b Data are from iKlvpin et al.l (j2002l ) . c The concentration of 
the NFW profile. d Data are from lHammer et al.1 ( 20071 ). e B/T is defined as the mass ratio of 
bulge to the total baryon mass. 



2.1. Dynamical model 



Following IBesla et al.l ( 120071 ) and iKallivayalil et al.l (120091 ). we constructed a model of 
the MW, M31 and of the LMC, the latter being considered as a point mass with a total mass 
of 2 x 10 10 M^ For both the MW and M31 we adopt a model consisting of a NFW halo 
(INavarro et al.lll997l ). a Hernquist bulge (IHernquistl Il990l ) and an exponential disk. Then 
the total gravitational potential of the galax y model is the sum-up of the three components 
jrlavashi et al.lhoOTl : Ishattow fc Loebll2009h : 



0b (r) 

M r ) 

<fa(r) 



4>b{r) + 4>d{r) + <ph{r), 
GM h 
r + a b ' 

GM d (l -e~^) 
r 

GM vir /r s 



ln(l + r/r s ) 



(1) 
(2) 

(3) 
(4) 



ln(l + c) — c/(l + c) r/r s 

where the subscript b, d, and h denote bulge, disk and halo, respectively; a b the scale length 
of Hernquist profile; Rd the scale length of disk; c and r s the concentration and scale length 
of the NFW profile. The parameters of the model are summarized in Table |5J The models 
are required to match the rotation curves of the MW and M31, respectively, see Fig. [TJ 



The equation of motion for each object can be written as: 



dt 2 



n 



d_ 

dfi 



(5) 



Then the trajectories of each object can be solved numerically using the standard method for 
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N-body simulation in barycentric frame. By choosing a small time step for the integration, 
it provides us an accuracy down to 0.1% over 10 Gyr, which is precise enough for our 
discussions. 

As mentioned in the intr oduction, nonsph erical halo of MW may have impacted to the 
trajectory of LMC as seen in (IBesla et al.ll2007l ). The cases of non-spherical MW halo can be 
studied by replacing r in 0h( r ) by r = \JR + z 2 /q 2 , where the cylindrical polar coordinates 



is ad opted, q characterizes the axis ratio of halo potential (IHayashi et al.l 120071 ; IBesla et al. 
20071 ). For q > 1 we refer to a prolate halo while q < 1 to an oblate halo. 



2.2. Results 

We uniformly sampled ~2000 possible M31 proper motions with the amplitude of 
[0,0.12] masyr -1 and the orientation [0,360] degrees on the sky. We define a reasonable 
solution by searching when the minimal distance between LMC and M31 can be less than 
50 kpc, enough close to be consistent with material ejected from an ancient merger, during 
the last 10 Gyrs. We do find a group of solution by using the LMC proper motion from 



Kallivayalil et al.l (|2006aJ). The solution of M31 proper motion for a spherical MW halo is: 

/iw = — 62 ± 18 /xasyr -1 , (6) 
/iN = — 25 ± 13 fiasyi^ 1 , 

where the error bar accounts for the error of LMC proper motions, and /iw, A*n are quoted in 
the equatorial system as usually used. It corresponds to f ra d=— 128km s -1 and f ta n=102km s _1 
for M31 relative to the MW. The averaged time since LMC was ejected from M31 is 
5.5 ± 1.4 Gyrs ago. Table 3 summarises the results after assuming different values for the 
axis ratio of the Milky Way potential. In this Table v ta , n , T tTave \ and v 50 are averaged values 
for all the trajectories that put the LMC at 50 kpc from M31 at lookback times indicated 
by Tt rav ei. At such a distance from M31 and for all solutions, the relative velocity of LMC 
to M31 was slightly higher than the escape velocity which is 408 kms -1 , consistently with 
expectations for material ejected from M31. Note that we define the escape velocity when 
an object arrives to the intergalactic space, i.e., 600 kpc from M31, between the Milky Way 
and M31. 

In Fig. H] we show the trajectories of M31 and LMC for the mean solution (q = 1). 
In the right panel of Fig. [2] we show the solutions of M31 proper motion varying with q 
that correspond to the average of all successful solutions. We have scanned a region of 
q = [0.5, 1.5] with span of 0.2, i.e., from oblate to prolate (see also Table |3]). Note that q 
is used in potential space, therefore the halo shape in density space would be even more 
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Table 3: Possible Solutions 



n a 


Mw 




^rad 


„. b,c 
^tan 


71 i 

- 1 travel 


^50 




( /"as yr" 


- 1 ) ( 


kms -1 ) 


( kms -1 ) 


(Gyr) 


( kms -1 ) 


0.5 


-83 ± 10 


01 ±08 


-127 


194 ± 44 


-4.3 ±0.5 


432 


0.7 


-77 ±14 - 


-10 ± 10 


-128 


160 ± 56 


-4.5 ±0.6 


428 


0.9 


-68 ±16 - 


-21 ± 14 


-128 


124 ± 53 


-5.1 ± 1.1 


423 


1.0 


-62 ±18 - 


-25 ±13 


-128 


106 ± 63 


-5.5 ± 1.4 


421 


1.1 


-56 ±20 - 


-28 ± 11 


-129 


89 ±68 


-6.2 ± 1.9 


418 


1.3 


-56 ±16 - 


-30 ± 09 


-129 


89 ±48 


-7.8 ±2.2 


417 


1.5 


-52 ±15 - 


-32 ± 09 


-129 


80 ±46 


-7.2 ±2.3 


417 



a The shape parameter of MW halo. b The velocities are given relative to the Milky Way. 
c The error bars actually delineate the solutions regions. d The velocity of LMC at 50 kpc 
to the M31 center. 

201 
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Fig. 2. — Left and Middle panels: Trajectories of LMC (red) and M31 (black) for the mean 
solution of a spherical MW halo (Eq. [6]). The open boxes indicate the time when LMC 
is closest to M31. Right panel: Possible M31 proper motions that satisfy our hypothesis. 
The gray dots are all the possible solutions, including by varying of MW shape, i.e., q, and 
accounting for the measurement errors of the LMC proper motion. The color dots connected 
by a line indicate the mean value for different q. The bars attached to each point indicate 
the solution region for each q. 



extreme for both prolate and oblate ones. 
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Discussion and Conclusion 



Could the Magellanic Clouds be ejected tidal dwarves from a previous major merger 
occurring at the M31 location? In this letter we simply demonstrate that it could be the 
case, within a reasonable range of parameters for both the Milky Way and M31. In our study 
we did not consider the effect of the Small Magellanic Cloud (SMC) to LMC as the former 
has a much smaller mass than the latter, and as such, our study of the LMC motion can 
realistically apply to the motion of both Magellanic Clouds together. We have also considered 
the possible impact of the gravitational pot ential of both the Virgo cluster and the Great 



Attractor (see e.g., iLynden-Bell et al.lll988l ). which are modelled by giant halos of NFW 



profile and find that they affect very marginally the LMC trajectories. The model presented 
here is possibly static assuming that only the M31 mass could have increase during the 
last 8 Gyr, through a major merger, thus neglecting some possible mass accretion of minor 
mergers. Assuming a smaller mass of M31 in the past (2 times less massive) would generate 
a larger travel time for the LMC (8 Gyr instead of 5.5 Gyr for our q = model). 

Besides this, we show that such an hypothesis is possible, although it does not demon- 
strate that it is indeed the case. To go beyond requires an estimate of several quanti- 
ties, especially the tangential velocity of M 31 and the axis ratio of the Milky Way poten- 
tial. Ivan der Marel fc Guhathakurtal ( 120081 ) have estimated tangential velocity of M31 to be 



/iw = —22 ±12 yuasyr -1 and /in = —11 ±10 /zasyr -1 , i.e. in the same direction than assumed 
in Eq. [6]), but with smaller amplitude. These estimates were based on the assumption that 
M31 satellites follow the motion of M31 through space. We have tested a null hypothesis for 
the tangential velocit y of M31 and find that it implie s a ve ry prolate Milky Way halo with 
q — 1.7. However the Ivan der Marel fc Guhathakurtal (120081 ) assumption may not hold in the 
case of an major merger in the past history of M31 because the orbital motions of its satel- 
l ite sy stem could be much more chaotic than expected by Ivan der Marel fc Guhathakurta 
( 120081 ). On the other hand, assuming a spherical halo for the Milky Way leads to values for 
the M31 proper motion which are larger than what is typically quoted, often on the basis 
of the timing argument. Possibly, the timing argument has to be re- formulated in a scheme 
for whic h there were more th an 2 bodies with a mass similar to the Milky Way, 6 Gyr ago 
(see e.g 



Hammer et al.ll2010f ). 



It is wiser to test our hypothesis by considering observable parameters that are not 
assuming a specific history for the Local Group satellite system. The result is somewhat 
troubling. First, by tracing back the LMC motion to M31, it is found that its relative 
velocity to M31 at 50 kpc is slightly above the escape velocity which is quite expected for 
tidal material ejected from a merger. Second, the travel time to reach the Milky Way is 
ranging from 4 to 8 Gyrs, depending on the axis ratio of the Milky Way potential (see Table 



-9 - 



3). lHammer et al.l (120071 ) estimated that if M31 have experienced a gaseous rich major 
merger, it should have occurred 5-8 Gyrs ago on the basis of the age of the M31 disk stars. 
Indeed in such an event most stars in the rebuilt disk should have ages slightly smaller than 
the merger look-back time. 

The origin of the Magellanic Clouds is still an enigma as they are the only blue, gas rich 
irregular in the immediate outskirts of the Milky Way. Our proposition has the advantage of 
explaining them in a consistent way, as being originating from the most massive body in the 
Local Group that show evidences for a very rich merger history. Future measurements of the 
M31 transverse velocity (possibly with GAIA) may confirm or infirm its validity. Further 
modeling of both LMC and SMC could be done to verify whe ther their internal struc tures 
(e.g. the LMC bar) and their star formation history (see e.g., lHarris fc Zaritskyl 120091 ) can 
be reproduced. Important tests of our hypothesis may come from better estimates of the 
da rk matte r cont ent of the LMC (could it be a tidal dwarf if it has a total mass as adopted 
by iPeebled (120101 )?) and from verifying whether it is consistent with the numerous features 
found in the outskirts of the Milky Way. 
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